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Abstract 

Background: Measures of attributable risk are an integral part of epidemiological analyses, particularly when aimed 
at the planning and evaluation of public health interventions. However, the current definition of such measures does 
not consider any temporal relationships between exposure and risk. In this contribution, we propose extended 
definitions of attributable risk within the framework of distributed lag non-linear models, an approach recently 
proposed for modelling delayed associations in either linear or non-linear exposure-response associations. 

Methods: We classify versions of attributable number and fraction expressed using either a forward or backward 
perspective. The former specifies the future burden due to a given exposure event, while the latter summarizes the 
current burden due to the set of exposure events experienced in the past. In addition, we illustrate how the 
components related to sub-ranges of the exposure can be separated. 

Results: We apply these methods for estimating the mortality risk attributable to outdoor temperature in two cities, 
London and Rome, using time series data for the periods 1 993-2006 and 1 992-201 0, respectively. The analysis 
provides estimates of the overall mortality burden attributable to temperature, and then computes the components 
attributable to cold and heat and then mild and extreme temperatures. 

Conclusions: These extended definitions of attributable risk account for the additional temporal dimension which 
characterizes exposure-response associations, providing more appropriate attributable measures in the presence of 
dependencies characterized by potentially complex temporal patterns. 
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Background 

Epidemiological studies usually rely on effect summaries 
based on ratio measures, such as relative risk, odds ratio 
or rate ratio, with the choice depending on the specific 
study design [1]. Although these measures are ideal for 
summarizing the association of interest, they offer limited 
information on the actual impact of the exposure. This 
information is critical for the planning and evaluation of 
public health interventions, and it is better provided by 
relative excess measures such as the attributable fraction 
(AF), or absolute excess measures such as the attributable 
number (AN). Steenland and Armstrong offer a thorough 
overview on the topic [2]. Here we generally refer to these 
summaries as attributable risk measures. 

Problems in the definition of these measures may arise 
in the presence of delayed associations, occurring when 
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an exposure generates a risk lasting well beyond the expo- 
sure period. Researchers in different fields, during the last 
thirty years, have proposed approaches for modelling this 
type of association [3-6]. In time series analysis, a popu- 
lar approach is based on distributed lag models (DLMs) 
[7,8], generalized to distributed lag non-linear models 
(DLNMs) when including non-linear exposure-response 
associations [9,10]. Recently, the DLNM framework has 
been extended beyond time series data for modelling such 
dependencies, defined exposure-lag-response associations, 
in different study designs [11]. However, attributable risk 
measures have not been developed for the DLM and 
DLNM class, with the result that their current definitions 
do not take into account the additional temporal struc- 
ture in the exposure-response association. Previous work 
investigated the issue in time series analysis, although 
without producing a general approach [12,13]. 

In this contribution, we extend this research and 
attempt to formalize the definition of attributable risk 
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measures within the DLNM modelling framework. In 
particular, we illustrate how complex non-linear and 
temporal patterns can be accounted for in the com- 
putation of the attributable risk. Also, we show how 
attributable components related to different exposure 
ranges can be separated. We propose an example on 
the estimation of the health burden attributable to tem- 
perature with time series data, a problem which has 
received quite a lot of interest recently due to climate 
change predictions. However, the approach can be easily 
applied to other exposure-lag-response associations. The 
method is implemented in simple functions developed 
within the R software package and provided as additional 
files. 

Methods 

Attributable risk measures 

A general definition of the attributable fraction AF* and 
number AN X for a given exposure x can be provided by: 

AF* = 1 - exp (-&) , (la) 
AN* = n ■ KY X , (lb) 

with n as the total number of cases. The parameter f) x 
used in Eq. (la) represents the risk associated with the 
exposure, and it usually corresponds to the logarithm 
of a ratio measure such as relative risk, relative rate or 
odds ratio. It is generally obtained from regression mod- 
els while adjusting for potential confounders. The general 
definition of f3 x used here refers to the association with 
a specific exposure intensity x compared to a reference 
value xq. For linear exposure-response relationships, the 
association can also be reported as /3 ■ x, where in this 
case fi refer to a unit increase in x. For binary variables 
reporting presence/absence of the exposure, Eq. (la) sim- 
plifies to AF = (RR — 1) /RR, with RR as relative risk, as 
reported by Steenland and Armstrong [2]. We keep the 
more general definition of /3 X , which is easily applicable to 
non-linear exposure-response relationships, throughout 
the manuscript. 

The theoretical nature of these effect measures is based 
on a counter/actual, where the observed condition is com- 
pared with a reference state which never occurred. This 
state postulates that the same population is followed in an 
identical situation where only the exposure level changes 
to the reference value xq. Typically, such a reference is rep- 
resented by the absence of association, meaning Xq = 0 
and p xo = 0. However, different counterfactual conditions 
can be used, for example a lower exposure which can be 
determined by an intervention. In this case the quantity 
P x can be simply re-parameterized as fi* = fi x — j3 XQ , and 
Eq. (1) still applies. 



Eq. (la) can be extended to define the risk attributable 
to multiple exposures X\, . . . ,x p : 

(P \ P 
■) = i -ri( i - AF *') - 

(2) 

with AN*!,...,* obtained by substituting Eq. (2) in Eq. (lb) 
[2]. For the specific form of Eq. (2), it should be noted 
that AFfl,...^ < AF^i + . . . + AF^, i.e. the sum of 
the attributable risk measured for individual exposures is 
usually higher than their concurrent attributable risk. 

A review of the DLNM modelling framework 

The basic idea underpinning the development of DLNMs 
is that the risk at time t can be described as the 
weighted sum of effects cumulated from a series of expo- 
sures Xt-i 0 , ■ ■ ■ ,Xt-L experienced in the past over the 
lag period £ = £q,...,L, with £q and L corresponding 
to minimum and maximum lags, respectively. The risk 
can be described by the function f(x), determining the 
exposure-response, and the function w{£), specifying the 
lag-response, related to the weights given to exposures 
at different lags £. These functions are combined in a 
bi-dimensional exposure-lag-response function f-w(x,£). 
Algebraically, the risk is defined by a function s(x, t; rj), 
written in terms of parameters r\ as: 

s(x,t;r))=f f-w(x t -i,£)d£ 
Ji 0 
L 

« J2f-">(xt-i,l)=vfl, t ri ■ (3) 

t=io 

The function s(x, t) is computed as the approximate 
integral of the exposure-lag-response function over the lag 
dimension, representing the cumulated risk over the lag 
period. The parameterization in the final step of Eq. (3) 
is obtained through a cross-basis, involving a tensor prod- 
uct between the basis chosen for f(x) and w(£), generating 
the transformed variables w x>t linearly combined with the 
parameters r). Simpler DLMs are defined by Eq. (3) by 
assuming f(x) as linear. Algebraic details and additional 
information are provided elsewhere [11]. The cross-basis 
is specified with a reference value xq used later as a cen- 
tering point for the function f(x), which is used to define 
the counterfactual condition. 

The complex parameterization of exposure-lag-response 
associations provided by Eq. (3) can be more easily inter- 
preted by computing effect summaries from the original 
parameters i). Specifically, the bi-dimensional exposure- 
lag-response risk surface modelled through /• w(x, £) can 
be expressed by a grid of effect summaries j3 Xl e, each inter- 
preted as the association with an exposure* at lag £ versus 
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the reference value xq. For a given time t, the cross-basis 
parameterization in (3) can be re-expressed as: 



L 

t=lo 



(4) 



This overall cumulative association is composed of the 
sum of contributions j3 Xi i from exposures Xt-i Q , ■ ■ -,Xt-L 
experienced within the lag period. Algebraic definitions 
have been previously provided [11]. 

Forward and backward perspectives 

The term fi x / for each intensity* can be interpreted using 
two complementary perspectives, illustrated graphically 
in Figure 1. From a forward standpoint, looking from 
current exposure to future risks, the terms /3 X / are the 
contributions from the exposure Xt occurring at time t to 
the risk at times t+£o> ■ ■ ■> t+L, identified by green circles. 
From a backward standpoint, looking from current risk to 
past exposures, the terms p x ,t are the contributions to the 
risk at time t from exposures Xt-t 0 , ■ ■ - ,Xt-L experienced 
at t — £q, . . . , t—L, identified by yellow squares. The under- 
lying curve in Figure 1 depicting such associations is called 
the lag-response curve related to a given exposure inten- 
sity x. The sum of these contributions over the whole lag 
period can be interpreted as the overall cumulative risk. 



attributable fraction b-AF^ t and number b-AN Xit at time 
t are obtained by substituting Eq. (4) in Eq. (2): 



(5a) 
(5b) 




b-AN^ = b-AF*, t • n t 



with n t as the number of cases at time t. This structure is 
consistent with the configuration of the regression model 
usually applied to fit the data, where the risk at time t is 
associated with lagged exposures at times t — I. The def- 
inition of backward attributable risk requires an extended 
version of the counterfactual condition accounting for the 
additional lag dimension: b-AN^ f and b-AF^j are inter- 
preted as the number of cases and the related fraction at 
time t attributable to past exposures to x in the period 
t — Iq, . . .,t — L, compared to a constant exposure Xq 
throughout the same period. 

An alternative version can be obtained using a forward 
perspective. Among other possible definitions, forward 
attributable number f-AN^j and fraction f-AF^f can be 
defined as: 



f-AF*,, = 1 - exp - 



(6a) 



Attributable risk from DLNMs 

The effect summaries provided above can be used for 
defining attributable risk measures within the DLNM 
framework. The idea is to treat the associations with expo- 
sures at different lags as independent contributions to the 
risk. A neat definition can be developed using a backward 
perspective, assuming the risk at time t as attributable 
to a series of exposure events in the past. The backward 



i-AN xt = f-AF* f • V — 

' ^ L - i o + l 

l=to 



(6b) 



This alternative version has some advantages if com- 
pared to the backward definition. First, the counterfactual 
condition is simpler: f-AF^ f and f-AN x , t are interpreted as 
the fraction and number of future cases in the period t + 
£o, . . . , t+L attributable to the single exposure x occurring 
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Figure 1 Conceptual model for the interpretation of exposure-lag-response associations: forward (left panel) and backward (right panel) 
perspectives. 
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at time t, compared to Xq. Moreover, the overall cumula- 
tive risk fS Xt} i for a given exposure Xt in (6a) is available 
also when the bi-dimensional exposure-lag-response is 
reduced to uni-dimensional exposure-response relation- 
ship, a step often needed in multi-site studies [14]. In 
contrast, all the lag-specific contributions are needed to 
compute Px t -i,t m (5a) for the backward counterpart. 

However, the forward version also has an important 
limitation, related to the fact that the contributions are 
associated to risks measured at different times. The 
attributable number f-AN^f in (6b) is computed by aver- 
aging the total counts experienced in the next Iq,...,L 
times, thus only approximating the lag structure of risks. 
This approximation is likely to produce some bias, which 
is expected as an underestimation of the attributable num- 
ber if compared to the backward version. 

Separating attributable components 

The definitions provided in Eq. (5)-(6) can be extended 
to separate the attributable components related to spe- 
cific exposures or exposure ranges. This will be used later 
in the example to single out the contributions from cold 
and heat in temperature-health associations. Let's define a 
range r =[l,h] between low and high exposure limits / and 
h. The definition of forward attributable number f-AN£ t 
and fraction f-AF£ t limited to exposures within the range 
r is clear-cut, as they are either equal to the quantities 
reported in Eq. (6) if x e r or zero otherwise. Adopting a 
backward perspective, a similar definition of h-AY r xt can 
be obtained by modifying Eq. (5a) as: 



b-AF; f = 1 - exp - 1 (*t-e e r) /3 Xt _ a , (7) 



simply selecting the risk contributions from past expo- 
sures included in the range r. The related attributable 
number b-AN£ t is computed by substituting Eq. (7) into 
Eq. (5b). Attributable components referring to different 
ranges can be summed up, as all are defined using the 
same counterfactual condition of a constant exposure 
Xt = Xo for the whole lag period I = Iq, . . . , L. 

The forward version has the additional advantage that 
for two non-overlapping ranges r\ and T2 the sum of 
the components is equal to the overall attributable risk, 
namely f-AF£ f +r2 = f-AF^+f-AF^. In contrast, adopting 
a backward perspective b-AF^ +r2 < b-AF^ f + b-AF^, as 
the risks are simultaneously computed for the same time t 
in the like of Eq. (2). 

Total attributable risk 

The attributable risk measures provided above can be 
computed for each of the i = l,...,m observations in 



a data set. An estimate of the total attributable number 
AN tot and fraction AV to t is provided by: 

m 

AN tot = J2 AN ^> ( 8a ' 

m 

AVtot = AN tof / J2 n k . (8b) 
;=i 

The equations above can be applied either to forward or 
backward attributable risk and to separate components, 
simply substituting the related attributable numbers in 
Eq. (8a). 

Computing uncertainty intervals 

Analytical formulae for confidence intervals of attri- 
butable risk measures are not easily produced [15], and 
this also applies to the extended versions developed here. 
Although approximated estimators have been proposed 
[15,16], in this context the most straightforward approach 
is to rely on interval estimation obtained empirically 
through Monte Carlo simulations [17,18]. Basically, we 
take random samples of the original parameters rj 
of the cross-basis in Eq. (3) from the assumed multi- 
variate normal distribution with point estimate t) and 
(co)variance matrix V(rj) derived from the regression 
model. These samples )/® are used to compute p^^ for I = 
io,...,L and each intensity x, empirically reconstructing 
the distributions of the attributable measures defined in 
Eq. (5)-(8). The related 2.5 th and 97.5 th percentiles of such 
distributions are interpreted as 95% empirical confidence 
intervals (eCI). 

Results 

The methods illustrated in the previous section are 
applied to estimate the all-cause mortality risk attributable 
to temperature, using daily time series from two cities, 
London and Rome, in the periods 1993-2006 and 1992- 
2010 respectively. R scripts and data implementing the 
method and partly replicating the results are provided as 
Additional files 1, 2, 3, 4, 5 and 6. 

Modelling strategy 

We fitted a standard time series Poisson model allowing 
for overdispersion, controlling for seasonal and long term 
trends and day of the week, using a 10 df/year spline and 
indicator variables, respectively. Model selection is still an 
issue of current research within the DLNM framework, 
although simulation studies indicate a good performance 
of methods based on the Akaike Information criterion 
(AIC) [11]. Considering the illustrative purpose of the 
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example, we selected a-priori the cross-basis function in 
Eq. (3) for representing the association between mean 
daily temperature and mortality, basing our choice on pre- 
vious analyses. Specifically, the cross-basis is composed 
of a quadratic B-spline with two equally-spaced knots as 
the exposure-response function f(x), and a natural cubic 
B-spline with three equally-spaced knots in the log-scale 
as the lag-response function w(l) over lags 0-25. 

In the specific case of temperature where a null expo- 
sure condition cannot be defined, a reasonable choice is 
to center the cross-basis in Eq. (3) to the temperature of 
minimum risk, as suggested in previous publications [13]. 
This optimal temperature corresponds to 20°C and 21°C 
for London and Rome respectively, and it represents the 
reference point xq for the computation of the attributable 
risk measures. These are obtained for the whole temper- 
ature range, and then for cold and heat contributions by 
separating the associations with temperatures lower or 
higher than In addition, the attributable components 
are separated further in mild and extreme cold and heat 
by selecting as cut-off values the 1 st and 99 th percentiles 
of city-specific distributions, corresponding to 0.4°C and 
23.7°C in London and 2.6°C and 28.6°C in Rome. 

We derived empirical confidence intervals for backward 
total attributable numbers and fractions, computed over- 
all and for separated components, by simulating 5,000 
samples from the assumed distribution of i). 



Risk attributable to temperature 

The estimated associations between temperature and all- 
cause mortality in the two cities are illustrated in Figure 2. 
The left panels show the bi-dimensional exposure-lag- 
response surfaces, while the right panels display the over- 
all cumulative exposure-response curves, interpreted as 
the risk cumulated over the entire lag period of 0-25 days. 
The associations are represented in the RR scale, with the 
centering point and the cut-off values for defining extreme 
cold and heat displayed as dotted and dashed vertical lines 
respectively. 

Table 1 reports the estimates of the total backward 
and forward attributable fraction b-AF tof and f-AF toi , 
with 95%eCI. The backward version indicates that overall 
13.59% and 12.58% deaths are attributable to temper- 
ature in London and Rome within the study periods, 
respectively. As expected, the corresponding estimates 
computed forward are affected by a certain degree of neg- 
ative bias associated to the averaging of future mortality 
within the lag period, as described above. Nonetheless, the 
difference is not substantial in this case. 

Cold, heat and extreme components 

The total backward attributable risk is then separated into 
components due to cold and hot temperatures, defined as 
those below and above the optimal temperature, respec- 
tively. The estimates, computed using Eq. (7), are reported 



Exposure-lag-response surface Overall cumulative association and temperature distribution 

London 1993-2006 London 1993-2006 




Temperature (C) 



Exposure-lag-response surface Overall cumulative association and temperature distribution 

Rome 1992-2010 Rome 1992-2010 




Temperature (C) 

Figure 2 Association between temperature and all-cause mortality. Left panels: 3-D graphs of the exposure-lag-response risk surface. Right 
panels: overall cumulative exposure-response associations with temperature distributions. London 1 993-2006 and Rome 1 992-201 0. 
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Table 1 Total mortality fraction (%) attributable to temperature, computed backward (b-AF fot ) and forward (f-AF fot ), 
reported as overall, hot and cold components with 95% empirical confidence intervals (eCl) 



Deaths Overall Cold Hot 



London 




b-AFfot 


13.59(10.04-17.09) 


12.95 (9.32-16.38) 


0.66 (0.52-0.80) 


845,215 






f-AF rot 


13.41 (9.72-16.87) 


12.84 (9.38-16.33) 


0.57 (0.45-0.68) 


Rome 




b-AF to t 


12.58(9.30-15.64) 


10.84 (7.37-14.23) 


1.74(1.12-2.37) 


395,691 






f-AF to f 


12.27(8.94-15.41) 


10.72 (7.19-14.00) 


1.55 (0.95-2.13) 



London 1993-2006 and Rome 1992-2010. 



in Table 1. The comparison of the two contributions 
clearly indicates that cold is responsible for most of the 
mortality attributable to temperature, with b-A¥ to t equal 
to 12.95% and 10.84%, compared to 0.66% and 1.74% for 
heat, in the two cities. Estimates of forward attributable 
risk are very similar, and as expected their sum is equal 
to the overall burden, differently than for the backward 
version. 

The analysis is extended by further separating the 
attributable components into contributions from mild 
and extreme temperatures, with results summarized in 
Table 2. For cold temperatures, the contribution from 
extreme days accounts only for a minimal part of the 
mortality burden, while still 12.3896-12.48% and 10.27%- 
10.37% of deaths in London and Rome, respectively, are 
attributable to mild cold days. This result is expected 
when inspecting the right panels of Figure 2, with the def- 
inition of mild cold including the majority of days in the 
series for both cities. 

In contrast, the comparison between the two cities is 
rather different for the components attributable to mild 
and extreme hot temperatures. In spite of the stronger risk 
in London, the attributable fraction is similar for extreme 
heat and even higher in Rome for mild heat (1.32%- 
1.45% versus 0.25%-0.33%). This apparent contradiction 
is explained by the different temperature distribution, and 
in particular the percentile corresponding to the opti- 
mal temperature, corresponding to 93.6 th and 72.5 th in 
London and Rome. This result suggests the hypothesis 
that the population in Rome is more adapted to the range 
of temperatures corresponding to extreme hot if com- 
pared to London, where the population experienced only 
a few days of unusually high temperatures. 



The harvesting paradox 

Accounting for the additional lag dimension in exposure- 
lag-response associations involves further complexities in 
the interpretation of attributable risk measures. We now 
focus our attention to the association with hot temper- 
ature in Rome. The left panel of Figure 3 shows the 
estimated lag-response curves at various temperatures, 
computed versus the reference optimal temperatures of 
21°C. The graph indicates a strong risk in the first lags, 
then followed by a protective association at longer lags. 
This pattern is consistent with the harvesting hypothesis, 
which assumes that the initial risk is partly discounted by 
the depletion of the pool of susceptible after an extreme 
heat event [19]. Also, the extent of harvesting seems more 
pronounced for extreme temperatures. 

This phenomenon has interesting implications. An 
example is offered by the right panel of Figure 3, illus- 
trating the estimated daily deaths b-AN X; f and i-AN Xt t 
attributable to heat, computed backward and forward for 
the first summer in the time series for Rome, with related 
temperature trend. As expected, a substantial number of 
deaths are attributable to temperatures above the opti- 
mal value, represented by the horizontal dotted line, in 
the period mid-July to mid-August. The trend of for- 
ward attributable deaths f-AN^ closely follows the daily 
temperatures, consistently with the definition of num- 
ber of deaths attributable to the temperature in day t 
cumulated in the next L days. In contrast, the back- 
ward attributable number b-AN^ decreases to zero and 
even becomes negative in late summer days, although 
the overall cumulative exposure-response in Figure 2 
(bottom-right panel) does not show a RR below 1 for any 
temperature. 



Table 2 Total mortality fraction (%) attributable to temperature, computed backward (b-AF tot ) and forward (f-AF fot ), 
reported as components from mild and extreme hot and cold contributions with 95% empirical confidence intervals (eCl) 



Extreme cold Mild cold Mild hot Extreme hot 



London 


b-AF to f 


0.55 (0.45-0.64) 


12.48 (8.86-15.88) 


0.31 (0.23-0.38) 


0.36 (0.29-0.43) 


f-AF to t 


0.47 (0.40-0.53) 


12.38(8.98-15.78) 


0.29 (0.22-0.35) 


0.28 (0.23-0.33) 


Rome 


b-AFfot 


0.59 (0.47-0.70) 


10.37(6.88-13.63) 


1.45 (0.89-2.01) 


0.33 (0.25-0.40) 


f-AF to , 


0.47 (0.39-0.54) 


10.27(6.69-13.50) 


1.32 (0.75-1.85) 


0.25 (0.19-0.30) 



London 1993-2006 and Rome 1992-2010. 
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Lag-response associations at different temperatures Daily attributable deaths and temperature 




Figure 3 Lag structure and harvesting paradox. Left panel: lag-response associations between various temperatures and all-cause mortality. 
Rome 1 992-2010. Right panel: daily number of deaths attributable to heat, computed forward (green circles) and backward (yellow squares), and 
temperature trend. Rome July-Sept 1992. 



This paradox is explained by the counterfactual con- 
dition associated with the backward perspective. Specif- 
ically, each b-AN x>t compares the association with the 
observed temperatures in the past L days to a constant 
exposure xq. In the presence of harvesting, the observed 
population becomes 'healthier' than the counterfactual 
population after a series of heat days, due to the deple- 
tion of the susceptible pool. This explains the negative 
attributable numbers for specific combinations of lagged 
exposures. This fact emphasises that harvesting should 
not be interpreted as a true protective association at 
longer lags, but rather as an artefact due to a change in 
the underlying population following a stress, which affects 
the counterfactual condition. This issue is relevant when 
using backward attributable risk measures b-AN^f and 
b-AF^t to assess the contribution of specific days. How- 
ever, similarly to the net overall cumulative risk, the total 
attributable number b-AN tot and fraction b-AN tot , pro- 
duced by Eq. (8) and reported in Tables 1 and 2, account 
for the discount by summing the contributions over the 
whole series. 

Discussion and conclusions 

In this contribution we illustrate an extended definition 
of attributable risk measures based on the DLNM frame- 
work. Consistently with this class of models, such a def- 
inition accounts for the complex pattern of potentially 
non-linear and delayed associations described through 
exposure-lag- response associations. 

Two alternative definitions of attributable risk are pro- 
posed, assuming backward or forward perspectives. The 
former provides more consistent estimators which nat- 
urally arise from the structure of the regression model, 



where distributed lag terms at times t — I contributes to 
the risk at time t. The forward attributable measures, in 
contrast, are affected from a negative bias related to the 
averaging of future counts, which nonetheless is likely to 
be relatively low. On the other hand, the forward ver- 
sion is well suited for separating the risk in components 
attributable to different ranges, as their sum matches the 
overall risk. Furthermore the forward perspective, looking 
from current exposure to future risk, seems more appro- 
priate for quantifying the health burden due to specific 
exposure occurrences, as it is based on a more coher- 
ent counterfactual condition. Corrections have been pro- 
posed in previous works on risk attributable to multiple 
exposures [20-22], and can be applied to the backward 
version. 

Strictly speaking, the definition given in Eq. la is 
interpreted as the attributable fraction among the sub- 
population of exposed subjects. In the setting of time 
series analysis for environmental stressors, the whole pop- 
ulation is usually considered as exposed, and this defini- 
tion can be more generally interpreted as the population 
attributable fraction. If only a subset is instead exposed, 
Eq. (5)-(8) can be easily extended using the equations 
proposed by Steenland and Armstrong [2] for population 
attributable risk. 

Previous papers suggested approaches for producing 
attributable risk from distributed lag models when applied 
to heat-mortality associations. Baccini and colleagues 
applied DLMs and computed attributable risk measures, 
specifically addressing the issue of harvesting [12]. Honda 
and colleagues illustrated an analysis on the mortality 
burden due to heat using DLNMs [13]. However both 
approaches are limited, as the former assumes a linear 
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threshold form of the exposure-response, while the latter 
averages the non-linear risk across the whole temperature 
range. In this paper we offer a formal and more consistent 
definition of such attributable risk measures. 

An advantage of the proposed method is the provision 
of estimates for separate components of the attributable 
risk, associated with different exposure ranges. In the spe- 
cific case of temperature-health associations, this allows 
the separation of attributable risks from cold and heat, and 
further from mild and extreme temperatures. The esti- 
mates reported in the example highlights how the simple 
analysis of exposure-response curves can be misleading in 
the attribution of risk, and that most of the mortality in the 
two cities is in fact attributable to mild cold temperatures, 
in spite of the relatively low RR. 

The availability of attributable risk measures, comple- 
mentary to estimates of exposure-response associations, 
is essential for the identification and planning of public 
health interventions. Their extension to exposure-lag- 
response associations allows the computation of such 
measures from dependencies showing potentially com- 
plex non-linear and temporal patterns. 

Additional files 



Additional file 1 : R script implementing the function to compute 
attributable risk measures. 

Additional file 2: Documentation of the attrdl function. 

Additional file 3: R script for fitting the regression models in the 
example. 

Additional file 4: R script for computing attributable risk in the 
example. 

Additional file 5: R script for producing the graphs in the example. 

Additional file 6: File with data used in the example. 
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